Observation of a Long- Wavelength Hosing Modulation of a High-Intensity Laser Pulse 

in Underdense Plasma 
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We report the first experimental observation of a long-wavelength hosing modulation of a high- 
intensity laser pulse. Side-view images of the scattered optical radiation at the fundamental wave- 
length of the laser reveal a transverse oscillation of the laser pulse during its propagation through 
underdense plasma. The wavelength of the oscillation Ahosing depends on the background plasma 
density n c and scales as Ahosing ~ n^ • Comparisons with an analytical model and 2-dimensional 
particle-in-cell simulations reveal that this laser hosing can be induced by a spatio-temporal asym- 
metry of the intensity distribution in the laser focus which can be caused by a misalignment of the 
parabolic focussing mirror or of the diffraction gratings in the pulse compressor. 
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Since the first theoretical prediction [l| and experimen- 
tal observation [2|-|4[ of laser-accelerated electron pulses 
showing a narrow energy spread, the prospect of a table- 
top laser-based particle accelerator has come closer to re- 
ality. Such an electron source would allow for the gener- 
ation of THz-radiation [5| , a high-intensity, narrowband 
femtosecond source of visible or X-ray radiation [fj, |7|, 
and possibly a free-electron laser in a university-scale lab- 
oratory. However, for all these applications a high degree 
of stability and reproducibility of the electron pulse is 
mandatory. This is difficult to achieve due to the inher- 
ently non-linear processes involved in the generation of 
quasi-monoenergetic electrons. Here, the ponderomotivc 
potential of the high-intensity laser pulse excites a plasma 
wave. Feedback between laser and plasma wave induces 
longitudinal and transverse modulation of the laser pulse 
and non-linear steepening of the wave. When the wave 
eventually breaks plasma electrons are injected into the 
electric field associated with the wave (the "wake field" ) 
|8| which can accelerate them to relativistic energies. 

To improve the stability of the electron beam, preion- 
ized plasma channels [9, [lfj] or gas- filled capillaries 11 | 
have been used. Additionally, improving the laser con- 
trast [12] and using two counter-propagation pulses 13] 
has increased the energy and pointing stability of the 
electrons. However, asymmetries in the spatio-temporal 
profile of the focussed laser pulse or of the initial plasma 
distribution have been shown for pulses significantly 
longer than the plasma wavelength to significantly af- 
fect the propagation of the laser through the plasma 
and eventually influence the interaction 14H17I]. Fur- 



thermore, the propagation of two separated focal spots 
can be interlinked via the plasma acting as a non-linear 
medium initiating a braided motion of the two focal spots 



] 18] . Thus it is crucial to gain as much information on the 
interaction process and about influences of experimental 
parameters that can significantly affect the interaction 
and eventually the electron-bunch generation. 

In this Letter, we show for the first time that also a 
laser pulse with a length of the order of the plasma wave- 
length exhibiting a distinctly asymmetric far-field inten- 
sity distribution can trigger a hosing modulation leading 
to transverse oscillations of the pulse during its propaga- 
tion through the plasma. The wavelength of this hosing 
strongly depends on the background plasma density. 

The experiments were carried out using the 10- TW 
TkSapphire laser system ASTRA at the Rutherford Ap- 
pleton Laboratory. 45-fs, 350-mJ laser pulses at a central 
wavelength of Al = 800 nm were focussed into the lead- 
ing edge of a supersonic He-gas jet produced by a conical 
nozzle having an exit diameter of 2 mm. This gas jet ex- 
hibited an almost flat-top longitudinal density profile as 
characterized before the experiment. The density in the 
plateau region could be adjusted by changing the back- 
ing pressure. Assuming full ionization the correspond- 
ing plasma densities were between n c = 4 x 10 18 and 
4 x 10 19 cm -3 . The horizontally polarized laser pulses 
were focussed by an //16 off-axis parabolic mirror to a 
focal spot with a diameter of 30 /im, the intensity aver- 
aged within this area was 1.2 x 10 18 W/cm 2 . The flu- 
ence profile of the focussed laser pulses was strongly at- 
tenuated and imaged onto a CCD camera. During the 
full-power shots, the pointing and an equivalent fluence 
distribution of the far-field of the laser could be moni- 
tored by imaging the leakage of the laser through the last 
plane mirror in front of the parabola using an f/20— lens 
and a CCD. Looking in the direction of laser polarization 
the interaction region in the plasma was imaged using a 



high-quality //2— lens with a magnification of 10 onto 
a 12-bit CCD camera equipped with an interference fil- 
ter transmitting at (800 ± 10) nm. The spatial resolution 
was close to 2 /zm, the field of view in the interaction 
region was 540 /im x 400 /im. MeV-electrons generated 
during the interaction could be detected by an electron 
beam monitor consisting of a CCD camera looking at a 
fluorescent screen positioned 480 mm behind the target 
and shielded by 20 mm thick aluminium making it sensi- 
tive to electrons with E^- m > HMcV. Furthermore, the 
spectra of the electrons could be measured using a high- 
resolution magnetic spectrometer equipped with imaging 
plate detectors. The electron spectra exhibited quasi- 
monoenergetic features. 

Fig. Q] shows four images of the scattered light 
recorded at the fundamental wavelength by the side- 
view imaging system. All images show a field-of-view 




FIG. 1: Side-view images of the electromagnetic radiation 
scattered from the plasma at Al- The images have a size of 
500 /imx 100/im and have been stretched vertically by a factor 
of 2. The laser pulse shows a hosing-type transverse oscilla- 
tion, highlighted by the white lines. The hosing wavelength 
continuosly increases with decreasing plasma density n e . 



of 500 /imx 100 /im. Note that they have been stretched 



vertically by a factor of 2. The initial plasma densities in 
these four images were n c — 2.8 (a), 2.4 (b), 2.0 (c), and 
1.6 x 10 19 cm~ 3 (d), respectively. All images have been 
filtered differently to record similar dynamic ranges for 
different levels of the scattered light. 

The laser pulse (propagating left to right) was focused 
into the leading edge of the gas jet located close to the 
left boundary of each image. The red ellipse in Fig.Q]a) 
depicts the size and approximate position of the vacuum 
laser focus. As the laser power is high enough to trigger 
relativistic self-focussing (the critical power for a density 
of n c = 2.8 x 10 19 cm" 3 is 1.2 TW, well below our laser 
power), the transverse diameter of the laser pulse de- 
creased during its propagation through the plasma. Note 
that since the radiation is only scattered from the plasma 
electrons at the instantaneous position of the laser pulse 
in the plasma, each image represents a history of the 
instantaneous position of the laser as a function of prop- 
agation distance. Hence, the horizontal axis represents 
both distance and time of propagation. 

Furthermore, it is obvious in all four images that the 
laser pulse does not propagate on a straight line but car- 
ries out transverse hosing-type oscillations. While the 
amplitude of the oscillation of a few ^m remains approx- 
imately constant both during the propagation and for 
the different background densities, the wavelength of this 
hosing modulation significantly increases with decreasing 
plasma density from (110±25) /im for the highest density 
up to (310±50) /im for the lowest. This is summarized in 
Fig. [2] Note that due to the decreasing plasma density 
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FIG. 2: Dependence of the hosing wavelength Ahosing on the 
electron density n c for the experiment (black solid squares) 
and for the results from the numerical simulation (green open 
diamonds). These results are compared to predictions by an 
analytical model scaling like Ahosing ~ n^ in 3D (red line) 
and like Ahosing ~ rQ ' in 2D (green line). 



the intensity of the scattered light decreases, since the 
scattering is directly proportional to the number of free 
electrons. It had decreased below our detection threshold 



for even lower densities. 

The long-wavelength hosing modulation observed in 
our experiment is very likely to be caused by the focal 
spot of the laser pulse exhibiting a spatio-temporal asym- 
metry in its intensity distribution relative to the phase 
fronts of the pulse. This asymmetry can arise due to a 
variety of reasons in an chirped pulse amplification laser 
system, such as by a slight misalignment of the parabolic 
focussing mirror, a slight misalignment of the diffraction 
gratings used to compress the pulse after amplification, 
pointing fluctuations of the laser pulse itself induced by 
mechanical vibrations or air turbulences or a combina- 
tion of these effects. Due to a feedback mechanism be- 




FIG. 3: Energy-density profile of the laser far-field at two 
different filter levels, a) and b). Snapshot of the initial, asym- 
metric laser-intensity profile used in the simulations, c). 



tween the intensity distribution of the laser pulse and the 
plasma density mediated by the local refractive index of 
the plasma, rj(r) = (1 — n e (r)/n cr ) 1 / 2 , and the pondero- 
motive potential of the laser, the propagation of the laser 
pulse is strongly influenced by the plasma density distri- 
bution. Here, n cr = oj^eorrie/e 2 = An 2 eom c c 2 / e 2 A 2 is the 
critical density of the plasma for the wavelength Al and 
frequency wl of the laser light, Eq is the dielectric con- 
stant and e and m e are the electron charge and mass, re- 
spectively. The gradient of the plasma's refractive index 
initiated by the head of the pulse via its ponderomotive 
potential will be asymmetric with respect to the bulk of 
the pulse, and will therefore provide an effective focussing 
"force" . This exerts a transverse, restoring force on the 
bulk of the pulse, causing this part of the pulse to exe- 
cute transverse oscillations. Typical, slightly asymmetric 
laser foci measured during the experiment with different 
filter levels are shown in Fig.[3]a) and b). Note that both 
images show a time-integrated measurement of the inten- 
sity distribution, i.e. the laser fluence in the focal plane. 
Under conditions where no such focal spot asymmetry 
was present, no hosing modulation could be detected. 

To verify this hypothesis, we carried out 2-dimensional 
Particle-In-Cell (2D-PIC) simulations using the code 
Osiris 19] starting with an asymmetric laser pulse. The 
spatio-temporal pulse intensity profile used in our simu- 
lations, as shown in Fig. [3] c) was achieved by overlap- 
ping two laser pulses of equal intensity, both having a 
radius of 10 //m (1/e 2 in intensity) mimicking the focal- 



spot asymmetry in the experiment. The two parts of the 
pulse were separated both in longitudinal (i.e. in propaga- 
tion direction) and in transverse direction by 6Al- They 
were overlapped at focus, so that the phase fronts were 
flat. This resulted in a smooth pulse profile despite being 
constructed of two pulses. The size of each part of the 
focal spot in the simulation was chosen as wq = 10 /im, 
which is slightly smaller than in the experiment to ac- 
count for the under-estimated self-focussing in 2D sim- 
ulations. Accordingly, clq = 1.5 was chosen higher than 
in the experiment. Fig. S](a)-(c) show the numerical re- 
sults obtained for three different densities. It is obvious 
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FIG. 4: Numerically simulated propagation of the laser 
pulse in the plasma. The images show the energy-density 
distribution the scattered light at the fundamental wave- 
length in transverse direction for electron plasma densities of 



(a), 1.75 X10 1 



(b),andl.0xl0 1 



(c) and (d). The corresponding hosing wavelengths are also 
indicated. No hosing occurs for a symmetric laser pulse (d). 



that the propagation of the laser pulse shows a transverse 
hosing modulation, with amplitude and wavelength sim- 
ilar to the experiment as also shown in Fig. [3] Fig. |U(d) 
shows the result for similar parameters as in (c) but for a 
pulse exhibiting no spatio-temporal asymmetry. In this 
case no hosing occurs. 

A heuristic description of the interaction can be found 
by assuming that the bulk of the pulse resides in the 
quasi-static wake generated by the head of the pulse due 
to its ponderomotive force. The general approach is to 
assume the pulse is weakly relativistic, so that a density 
depression but not complete cavitation is generated by 
the pulse-front. The bulk of the pulse, displaced trans- 
versely to the propagation direction with respect to the 



head of the pulse, experiences the refractive index gradi- 
ent associated with the density depression. This acts to 
direct the pulse energy towards the axis. 

By considering the bending of phase fronts of the bulk 
of the pulse analogous to [20| , the position of the centroid 
of the bulk of the pulse, (r) , with respect to the radially 
varying plasma density induced by the front of the pulse, 
neglecting longitudinal dynamics, can be expressed as 



dt 2 



(r) 



~2dr 



6n c 
n cr 



(1) 



Assuming an approximate balance between the pon- 
deromotive and wakefield potentials for a weakly rel- 
ativistic pulse, the density perturbation can be ex- 
pressed as 5n e /n CI « (c 2 /w£)V 2 a 2 /4. If the held 
strength is distributed as a(r) 2 = ciq(wq/w 2 ) cos(nr/2w), 
where wq and clq = eEo/(m e cu>L) are the initial pulse- 
front radius and the normalized amplitude of the vec- 
tor potential, respectively (Eq is the amplitude of 
the laser's electric field), and w is the radius of the 
front of the pulse after self-focussing, then Sn c /n CI w 
~(aoW§7r 2 c 2 /16o; 2 j w 4 )cos(7rr/2w). The expression for 
the field strength neglects absorption of pulse energy 
which conserves laser power, a^uio = aw. We will as- 
sume that the pulse asymmetry lies in the x — z plane so 
that the oscillations will occur in the ^-direction. Hence, 
for the bulk of the pulse the equation for the centroid 
position is now 



dt 2 ' 
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To first order in (x), this is a harmonic oscillator equa- 
tion with a frequency closing ~ aoWoTT 2 c 2 /(8y/2uJLW 3 ). 
Using this frequency of oscillation, the expression yields 
Ahosing ~ 1 6 v / 2wlw 3 /(oo wo ttc). For a weakly relativistic 
pulse, the pulse self-focusses to a radius of approximately 
half the plasma wavelength, w « ttc/uj p [2lf , leading to 
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(3) 



For initial pulse radius and field strength Wq — 15 /im, 
a = 0.75, respectively, and Al = 0.8 //m, a hosing wave- 
length of Ahosing = 0.32^m- (n cr /n c ) 3 / 2 is obtained. This 
agrees well with the experiment as shown by the red 
line in Fig. [2] For the simulations, in a 2D slab geom- 
etry, the expression for constant power is cl^wq = a 2 w. 
This modifies the hosing wavelength to be Ahosing. 2D ~ 
8A L ' /(aoy / u>o)(n cr /n ) 5 / 4 , which explains the slightly 
different scaling. Using iuo=10/im, oo=l-5 gives the 
dashed green line in Fig. [2j which shows excellent agree- 
ment between the model and the 2D-simulation data. 

In conclusion, we have observed the transverse hosing- 
type oscillation of an ultra-short laser pulse during its 



propagation through underdense plasma. Analytical es- 
timates and 2D-PIC simulations have shown that such a 
transverse oscillation can be caused by a spatio-temporal 
asymmetry of the intensity of the laser pulse. For ex- 
periments towards the generation and application of col- 
limated, quasi-monoenergetic electron pulses where sta- 
bility is an important issue, such transverse oscillations 
of the laser pulse are likely to influence the electron accel- 
eration process. To improve the stability of the electron 
pulse it is necessary to carefully control and monitor the 
spatio-temporal characteristic of the driving laser pulse. 
It is conceivable that the hosing modulation will become 
more pronounced at higher laser powers as in this case 
the sensitivity to phase front tilts is likely to be higher. 
Deliberately exciting hosing of this form could also have 
a significant effect on x-rays produced by betatron oscil- 
lations which has been shown to be sensitive to focal spot 
asymmetries [22j. 
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